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ABSTRACT 

Magnetic diffusion in accretion flows changes the structure and angular momentum 
of the accreting material. We present two power law similarity solutions for flattened 
accretion flows in the presence of magnetic diffusion: a secularly-evolving Keplerian 
disc and a magnetically-diluted free fall onto the central object. The influence of Hall 
diffusion on the solutions is evident even when this is small compared to ambipolar 
and Ohmic diffusion, as the surface density, accretion rate and angular momentum in 
the flow all depend upon the product ^jj(B-fi), and the inclusion of Hall diffusion may 
be the solution to the magnetic braking catastrophe of star formation simulations. 
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1 INTRODUCTION 

Accretion flows onto protostars are flattened and braked by 
magnetic fields. The transportation of angular momentum 
from the flow to the envelope is facilitated by the bending 
of the magnetic field lines and inhibited by the drift of those 
lines against the flowQ. In magnetohydrodynamical (MHD) 
simulations of star formation it has been shown that it is 
possible to remove all of the angular momentum from the 
accretion flow using a weak field, with the consequence that 
no ac cretion disc forms (e.g. lPrice fc Bat e 2007; Mcllor Tfc Lil 
2008). Magnetic diffusio n reduces the braking so that a pro- 
tostellar disc forms (e.g. lKrasnopolskv fc Ko nigl 2002, here- 
after KK02; Mellon & Li 2009); these discs are precursors 
to those in which planet formation occurs. 

The amount of magnetic braking affecting an accretion 
flow depends on the coupling of the largely neutral medium 
to the magnetic field, which is facilitated by collisions be- 
tween the neutral and charged particles that transmit the 
Lorentz force to neutrals. In most simulations the magnetic 
field behaviour has been approximated by ideal MHD (where 
the magnetic field is frozen into the neutral particles), how- 
ever this approximation breaks down as the density in the 
disc increases. The relative drifts of different charged species 
with respect to the neutral particles delineate three conduc- 
tivity regimes. Ohrrnc (resistive) and ambipolar diffusion, 
which behave similarly in a thin disc (to a first order ap- 
proximation), have b een shown to reduce the effectiveness 
of magnetic braki ng l|Shu et all [20061 : iMachida et al.ll2007l : 
iMellon fc lH|20Q9| ). Hall diffusion occurs when the degree of 



coupling between the different charged species and the neu- 
trals varies, and is expected to dom i nate in large regions in 
accretion flows (|Sano fc Stondl2002l : rWard"ie! r2007). For cer- 
tain field configurations Hall diffusion can in principle cause 
the magnetic field to be accreted faster than th e neutral fluid 
jWardle fc Nelll999l : [Braiding fc Wardle!l2012l ). 

The role of Hall diffusion in star formation and accre- 
tion discs is neither fully understood nor explored due to the 
difficulty of performing n umerical simulations (although this 
is sta rting to change, e.g. iLi et al.ll201lj : iKrasnopolskv et al.l 
2011), highlighting the need for simpler analytic models that 
demonstrate the importance of the Hall effect. In this paper 
we present two power law similarity solutions to the MHD 
equations for an isothermal thin disc with Hall, ambipo- 
lar and Ohmic diffusion: one in which the disc is Keplerian 
and a second in which the collapsing material spirals onto 
the central object without disc formation. We discuss the 
implications and the physics controlling these remarkable 
solutions. 



2 FORMULATION 

The magnetohydrodynamic equations for an isothermal sys- 
tem are given by 



! + v.( P v) = o, 



(1) 



p^+p(V-V)V = -VP-pV$ + JxB, (2) 
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This is only strictly true for ambipolar and Ohmic diffusion. 



V 2 $ = AirGp, 



(3) 
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V • B = 0, 



(4) 



and 
dB 



V x (V x B) 



-Vx [j 7 (VxB)+7 7 H(VxB)xB + r ? A(VxB) ± ], (5) 

where p is the gas density, V the velocity field, P the mid- 
plane gas pressure, <3> the gravitational potential (defined as 
g = — V$ where g is the gravitational field), G the grav- 
itational constant, B the magnetic field and J the current 
density defined as J = c(V x B)/47r where c is the speed 
of light. The diffusion coefficients for the Ohmic, Hall and 
ambipolar terms in the induction equation are t], t)h and tja 
respectively. 

We adopt cylindrical coordinates and assume that the 
disc is axisymmetric and thin. The velocity field, V(r), has 
both radial and azimuthal field components, the latter giv- 
ing rise to the specific angular momentum J = rV^, as the 
material is assumed to settle rapidly to the disc midplane. 
We assume that the disc is threaded by an open magnetic 
field that is symmetric about the midplane, defining B z as 
the vertical field component at the midplane and B rs and 
B^s as the components at the surface z = H, where H(r) 
is the half-thickness of the disc. We neglect any mass loss 
due to a disc wind. As the disc is isothermal, the pressure at 
the midplane is given by P = pc 2 , where c s is the constant 
isothermal sound speed, typically taken to be 0.19 km s~ . 

Equations HH 5 ] are vertically-averaged as in 
iBraiding fc Wardle] (|2012l , following KK02) by assum- 
ing the disc is thin and integrating over z, which allows 
us to discard any terms of order H/r. The density, radial 
velocity, azimuthal velocity and radial gravity are approxi- 
mated as being constant with height, as are the quantities 
77, rjH I B and t/a/B 2 . The radial and azimuthal magnetic 
field components are taken to scale linearly with height, so 
that their values at the disc surface (B rs and B^ s ) may be 
used to characterise the field. 

The equations are further simplified by employing 
monopole expressions for B rs and g r : 

*(r,f) 



B rs — 



and g r = 



2-KT 2 

GM(r, t) 



(6) 
(7) 



where M and \& are the mass and magnetic flux enclosed 
within a radius r. These describe the disc sufficiently well 
that a more complicated iterativ e method of calculating t he 
field and gravity is unnecessary (|Contopoulos et al.| [l998). 

The azimuthal field component is calculated by bal- 
ancing the torques on the disc from rotation and magnetic 
braking, which transports angular momentum from the disc 
to the lo w-density external medium into which the field 
is frozen ( Bas ufe Mouschoviasl Il994r i . Assuming that the 
background angular frequency is small and that the exter- 
nal Alfven wave speed is constant and parameterised by 
VA.ext = Cs/ct, the azimuthal field component is then 



B, , 



VhB,. 



HB 



*a T) A B z 



2 c a B 2 H 



;5B Z 



(8) 



jBraiding fc WardkH20l3 l The field is capped as a way of 
representing the many magnetohydrodynamical instabilities 
and processes that might prevent the azimuthal field com- 
ponent from exceeding the poloidal component such as in- 
ternal kinks, turbulence and the magnetorotational instabil- 
ity (MRI). The cap 8=1 corresponds to the typical value 
of \B l f >s \/B z obtained when the vertical angular momentum 
transport is dominated by a centrifugally-driven disc wind 
(KK02). 

The vertically-averaged equations are then: 



dV r , T/ 9V r 



c|gE B z B rs J 2 
E dr 2ttE r 3 



2£ + V — = rB * B 4>° 
dt r dr 2ttE : 



(9) 



(10) 



(11) 



Ec s = n GE2 + GM^H + 1 {Bls + B%) (12) 



2H 2 



4r 3 8tt 



and 



H 9 * _ r uv n r,R rrtH r r rr]A n r 2 . 
2^~Qj~ —rnV r B z — rjB rs B —B z B rt>s -^B va B z \ 

(13) 

E is the column density, defined as E = 2Hp, and M c is the 
mass of the central s tar. These equation s are a simplified 
set of those used in IBraiding fc Wardle! l|2012T ). removing 
the HdB z /dr terms that were used to refine the structure 
of the disc and shocks in the full similarity solutions. 

The Ohmic and ambipolar diffusion terms scale to- 
gether, to a zeroth-order approximation, as they possess 
a similar dependence upon B and appear in the induction 
equation multiplied with the same field component. Ohmic 
diffusion is expected to be more important in the inner re- 
gions where the density is high , while ambip olar diffusion 
dominates in the outer regions (| Wardle! l2007f ). As the field 
within the disc is effectively vertical, ambipolar and Ohmic 
diffusion influence the field drift in the same direction, and 
only one term is required to study the change in disc be- 
haviour introduced by Hall diffusion. We then describe am- 
bipolar and Ohmic diffusion using the Pedersen diffusivity, 



VP = VA + V, 



(14) 



which we parameterise by the constant nondimensional Ped- 
ersen diffusion parameter, fjp (described in Appendix A). 

Similarly, we define a nondimensional Hall diffusion pa- 
rameter f]H such that the Hall coefficient r]H scales with the 
surface density and thickness of the thin disc in a similar 
way to the Pedersen coefficient (see Appendix A). The direc- 
tion of Hall diffusion depends upon the product r;^(B • Q), 
wher e the sign of t)h de pends on microphysics within the 
disc l|Wardle fc Ng| 11999); in our simulations we vary the 
sign of f\H in order to examine the effect of reversing the 
magnetic field. By studying the range of parameters that 
give stable accretion flow solutions it may be possible to 
place constraints on their physical values in observed accre- 
tion systems. 

We look for solutions to the equations that take the form 
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of power laws with respect to a similarity variable x — r/c s t 
in the limit that x — > 0. These may then be thought of as 
valid in the limits r — > or t — > oo, that is, the innermost 
regions or late stages of the accretion flow. Only two phys- 
ical solutions exist to these equations (see Appendix [XJ: a 
secularly-evolving Keplerian disc and a magnetically-diluted 
near-free fall collapse. We present these in the following two 
sections. 



3 KEPLERIAN DISC SOLUTION 

The first solution is a Keplerian disc, where the material 
is supported against gravity by rotation (V^> = ±^/GM /r; 
the z-axis is defined such that B z is always positive) and 
accretion onto the central mass is slow. The disc is described 
by the simplified set of equations: 



J 

d J _ rB$ s B z 
£fr ~~ 2ttT,V t ' 

HV r + ^-B^s H — ^-B ra = 0, 
B B 

<;Mc -H 2 +irGZH-c 2 s =0, 



2r 3 



B r 
Bo 



M = constant, 
J = rV^, 
4 . 



-B z 



SB Z 



8 2 
and $ = — ivr B z \ 



(15) 
(16) 
(17) 

(18) 

(19) 
(20) 

(21) 

(22) 

(23) 



the other terms in Equations [51 4131 particularly the deriva- 
tives with respect to time, are negligible in this limit. The 
induction equation (|17[1 can be written quite simply as 
V r + Vbt = 0, that is, the inward radial velocity of the fluid 
is balanced by the drift of the magnetic field with respect 
to the gas (Vbt)- Any accretion through the centrifugally- 
supported disc is regulated by the outward diffusion of the 
magnetic field against the flow; the magnetic diffusion then 
determines the column density E. The accretion rate M and 
the direction of rotation in the disc are determined by the 
magnetic torque B^ S B Z = —SB 2 , and are not affected by a 
global reversal of the magnetic field. 

Hall diffusion can either add to or act against the Ped- 
ersen diffusion in the disc, as r\n can be either positive or 
negative, while r\ and t/a are both positive. The cap on the 
azimuthal field, 5, can have either sign, creating counter ro- 
tating solutions when 8 is negative. Hall diffusion acts in 
the opposite direction in the counter rotating solutions (as 
the Hall term in Equation[T7]has the opposite sign), causing 
there to be symmetry between the clockwise solutions with 
t\h positive and large and the anticlockwise solutions with 
large negative r\H- 

The power law solution is given by: 



M 



~~G~ 

m 
cri 



1/2 



M: 
E : 



mc a 

5/2 ,1/2 

2-kG rW'- 



(24) 
(25) 



B z 
H 



= ± 



rnc-H 



(26) 



r 3/2 

hi- 



„, ftr(mcg) 3 /*_3/ V /* (27) 



3^/2\S\G 



>c a t 



(28) 



where the nondimensional constants characterising the col- 
lapse are m, the nondimensional infall rate onto the central 
star, 5, the artificial cap placed on the azimuthal magnetic 
field component, and fju and f/p, the Hall and Pedersen dif- 
fusion coefficients. The sign of J and V$ matches that of S, 
as the azimuthal drift of the fluid against rotation causes the 
field lines to bend at the disc surface. The magnitude of the 
cap on the magnetic braking, not its direction, influences 
the vertical field and the surface density. Finally, o\ and h\ 
are constants, determined by the expressions 



hi 

and <Ji — 



m[l + (f/28)*] 
V2?nf 



(29) 



(30) 



2|5| v /(2«5//) 2 + l 
(from Equations ! 17l and ll8p . where / is given by the equation 



f=±fj P -6fj a J^ 



S 2 . 



(31) 



This coefficient must be positive for accretion to occur 
through the disc — the opposite sign corresponds to a solu- 
tion in which the field is diffusing inwards against the out- 
wards flow of matter. For the typical value of 5 = 1, the 
condition / ^ 1 becomes 



fjp ^ 



17 

8 ' 



(32) 



creating a forbidden region of parameter space, shaded in 
Fig. [T] and bound by the solid line / = E = 0. 

Fig-HJiUustrates the two-dimensional region of magnetic 
diffusion parameter space that gives feasible values of the 
surface density E at r = 1 au when t — 10 4 years, c s = 0.19 
km s _1 and 5 = 1 for solutions with an accretion rate of 
M c = 10" 5 M yr" 1 (corresponding to B z = 1.15 G). The 
dotted line represents a surface density at 1 au of E = 1700 g 
cm -2 ; this is the value from t he minimum mass solar nebula 
model l|Weidenschillinelll977h in which the surface density of 
the solar nebula is estimated by adding sufficient hydrogen 
and helium to the solid bodies in the solar system to recover 
standard interstellar abundances, and spreading this mate- 
rial smoothly into a disc. The dashed lines correspond to 
higher surface densities that are more like those expected to 
occur in protostellar discs. 

The radial scaling of the surface density (E oc r~ z ^ 2 ) in 
this solu tion is that expected from the minimum mass solar 
nebular (|Weidenschillinell 19771) and other simulations of pro - 
tostellar discs (e.g. lTsuribelll999l : IVorobvov fc Basu] [2009): 
and the magnetic field scaling also matches that from theory, 
particularly of discs that support disc winds (B z oc r -5 ^ 4 ; 
IBlandford fc Pavne|[l982l ). The azimuthal field component 
blows up with respect to the other field components in this 
small x limit because the azimuthal magnetic field drift 
speed is slow compared to the Keplerian speed. The model 
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Figure 1. Exploration of the relationship between the Hall and 
Pedersen diffusivities r\n and fjp in the Keplerian disc solution for 
various disc surface densities £ at r = 1 au and t = 10 4 yr, when 
the azimuthal field cap 5 = 1, the sound speed c B = 0.19 km s _1 
and the accretion rate is M = 10 -5 Mq yr~ 1 (corresponding to 
B z = 1.15 G). The solid black line is £ = 0, and the shaded region 
beneath has no solutions describing an accreting disc rotating 
in the clockwise direction. The dashed lines correspond to £ = 
1700, 3000 and 10000 g cm -2 respectively; the equivalent counter 
rotating solutions are reflected about fjjj = 0. 



adopted for the vertical angular momentum transport is un- 
able to properly account for the effects of magnetic braking 
in the small x limit, so B^ s must be capped. If a different 
scaling for t\h were adopted then Hall diffusion could act 
to limit Btfts such that the cap becomes unnecessary; this 
should be examined in future disc studies. 

As is clear from Fig. [1] there is no solution describing 
a disc with purely Hall diffusion where the Hall diffusion 
parameter is positive, as the positive Hall parameter acts to 
restrain the effects of ambipolar and Ohmic diffusion when 
the disc is rotating in the clockwise direction. However, when 
the Hall diffusion parameter is negative (corresponding to a 
reversal of the magnetic field) it acts in the same radial direc- 
tion as the Pedersen diffusion, enhancing the radial diffusion 
of the magnetic field. When the disc is counter rotating, the 
opposite holds true, so that the solutions in Fig. [1] are re- 
flected across the vertical line fjn = 0. In the case of pure 
ambipolar diffusion the field moves inward slower than the 
neutral particles and the solution reduces to the fiducial so- 
lution of KK02 (their equations 51-57). 

Fig.[2]plots E against B z and M, showing the influence 
of Hall diffusion on E at 1 au of a disc with fjp = 8 = 1. The 
surface density increases with the negative Hall parameter 
to create heavier discs, while a positive value of t)h causes E 
to drop off dramatically as t)h approaches the limit of / ^ 0. 
Once more the opposite behaviour holds true when the disc 
is counter rotating. 
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Figure 2. Disc surface density £ against the vertical magnetic 
field component B z and central mass accretion rate M for dif- 
ferent values of the Hall diffusion parameter f\n for a Keplerian 
disc rotating in the clockwise direction with Pedersen parameter 
f] P = S = 1.0 and c = 0.19 km s _1 , at r = 1 au and t = 10000 yr. 
The solid line corresponds to f\n = 0; the dashed lines above are 
f]pi = —0.2, —0.5, —1.0 and —2 respectively; and the dotted lower 
lines are fju = 0.2, 0.4 and 0.6. 



Equations ll5II20l mav be written in a more physical form 
by writing the variables as functions of B z , r, t, the diffusion 
parameter / and the azimuthal field cap. The magnetic field 
components are then given by Equations I21H23I while the 
other variables are: 



M = 
M = 



f(26f 



Bt 



V Gt 



Bt 



V G 



1/3 



1/3 



Csf 

27r|5| 



M 



2G[(2<5//) 2 + 1] 



-3/2 



H = 



J 

and V r 



G[l + (f/2Sr]M 
(28GB 2 z r 4 t) 1/3 , 
(28GB 2 z rt) 1/3 
SBlr 



3/2 



ttV^E 



(33) 
(34) 
(35) 

(36) 

(37) 
(38) 

(39) 



This model description shows more clearly that the direction 
of rotation in the disc is tied to the direction of magnetic 
braking and emphasizes the build up of the central mass 
occurs through accretion against magnetic field diffusion. 

Alternatively, the disc variables can be thought of as 
functions of M and M, by treating t as the ratio of M to M 
rather than the age of the system. The solution can then be 
regarded as a steady state Keplerian disc, where the vari- 
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ables are defined by 



B z = 



/ GMM i 
V (25) 2 r 5 



J = VGMr 



V 4> = V K = 



GM 



and V r = 



M 



(40) 
(41) 
(42) 

fcrEr ' (43) 
E and H are unchanged from Equations 1351 and 1361 and the 
field components are as in Equations 12 lt[23l 

The amount of Hall diffusion in the flow helps deter- 
mine how much matter is accreted onto the central object 
compared with that remaining in the disc; decreasing the 
Hall parameter f/n in comparison to a constant i)p causes 
the surface density become large as the radial field diffusion 
slows the rate of infall. The ambipolar diffusion-only fiducial 
solution of KK02 has an accretion rate of M c = 7.66 x 1CP 6 
Mq yr _1 , and the surface density at 1 au is 1310 g cm -2 . 
Adding a positive value of t)h to this solution causes the 
surface density to decrease, which is problematic if one ex- 
pects to form large planets in a protostellar disc, however, 
if the Hall parameter were negative then the surface den- 
sity would be raised to a more realistic value. This again is 
only true for solutions rotating in the clockwise direction; 
however, as the sign of t]h depends on microphysics in the 
disc, the general point that planetary system architecture 
depends on the product r]n(B ■ fi) stands. 



4 FREE FALL SOLUTION 

The second similarity solution describes the infall when the 
magnetic braking is very efficient at removing angular mo- 
mentum from the flow. In this case there is very little angular 
momentum remaining and the reduced centrifugal support 
inhibits disc formation, so that the collapsing flow becomes 
a supersonic magnetically-diluted free fall onto the central 
protostellar mass. This solution is representative of the mag- 
netic braking catastrophe that affects many numerical sim- 
ulations of gravitational collapse. The collapse is described 
by the equations: 

dV r 

= 9r, 



dr 
dJ_ 

dr 



rB 4>a B z 
2tt£V; 



V r H + ^—B^s + "%p;B r 

t> D 



B Z = 0, 



GM a (-B; 

2r 3 4tt£ 
M — constant 
J = rV^, 
B rs = B, 



H-ct = 0, 



* = 2nrB z 



and Bd 



—f]H 
VP 



B ; 5B Z 



(44) 
(45) 
(46) 

(47) 

(48) 
(49) 
(50) 
(51) 

(52) 



as in the Keplerian disc solution the induction equation takes 
the simplified form V r + Vp.r = 0, however here the magnetic 
diffusion is assumed to be large and dynamically important. 



In this similarity solution any rotation remaining in the flow 
is that induced by the magnetic "braking" , which can cause 
rotation by Hall diffusion of the field lines tied to the elec- 
trons in the azimuthal direction, which creates a rotational 
torque on the neutrals and grains as they fall inward rapidly. 

The dimensional form of this similarity solution is given 
by the complete set of fluid variables: 



mc a 

IT 



M 

Vr 

V4, = jl,2C 

B 



t, 



2rhcit 



2 u 
C s 0zl,2 



Br 



and 



B z 



Bd>s = — min 



M : 

E : 
J 

* = 
H = 



~~G 

1 

~Tg 



1 mc% 

: ji,2C s r, 
2%c 2 a b zl)2 

Ql/2 r > 
r 3/2 

hl,2 ■ 



-rjH 



B ; 8B Z 



VP 



(53) 

(54) 
(55) 
(56) 

(57) 



(58) 



The sign of the Hall diffusion coefficient and the cap on B^s 
must be taken into account when calculating the minimum 
of Equation 1581 as this should be an absolute minimum. We 
define the z direction such that B z is always positive, and as 
the radial velocity is negative then conservation of angular 
momentum (Equation I49p requires B$ s and J to have the 
opposite sign. Any twisting of the field by Hall diffusion is 
balanced by spinning the fluid in the opposite direction. 

Given that there are two possible values for B$ a , de- 
pending upon the chosen values of the parameters v~h,p and 
S, there are two possible coefficients for B z , J and H . These 
are referenced by the subscripts 1 and 2 depending on which 
term in Equation [58] is smaller. The first set of coefficients 
occur when 



ffp-Vs 



Bs a — 



so that the azimuthal field is given by 

y/2y H B z 

The vertical field coefficient b z \ is then the single 
of the polynomial 



8 m , 6 m 

Oxl - n ~n c °zl - 



2vlfi 



= 0, 



where /1 is a constant defined by the magnetic 
parameters: 

, V 2 p + V 2 h 

,fl = ^2"' 

Tp - Vh 

The coefficients for j and h are given by 



and 



hi = 



Ji = 



fibli 



™ \fW P - vl 



-1 + 4/1- 



4m 2 

Wii 



(59) 

(60) 
real root 

(61) 
diffusion 

(62) 
(63) 
(64) 
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positive real root of the equation 



Figure 3. The relationship between the Hall and Pedersen dif- 
fusion parameters t}h and fjp in the free fall solution, for various 
midplane magnetic field strengths B z at r = 1 au and t = 10 4 
yr, where the azimuthal field cap 5=1 and the sound speed 
c a = 0.19 km s _1 . We take the central mass accretion rate to be 
M = 10 — 5 Mq yr -1 (corresponding to a surface density at r = 1 
au of E = 5.03 g cm -2 ). The solid curves correspond to negative 
rotation, the dashed ones positive, and the boundary between the 
two values of B<p s is the dot-dashed line rjp = fiH\/2 + S 2 /S. 



When there is only ambipolar diffusion, this solution reduces 
to the nondimensional asymptotic behaviour of the "strong 
braking" collapse of KK02 (their equations 66-71). In their 
solution fi = 1 and ji « 0, so that the angular momen- 
tum and the azimuthal field component are reduced to a 
small plateau value as the strong magnetic braking removes 
almost all of the angular momentum early in the collapse 
(see their fig. 10). Their similarity solutions tended towards 
this behaviour when the magnetic braking parameter a was 
large, even though it does not appear in the coefficients of 
the power law solutions. 

This set of coefficients corresponds to the curves in Fig. 
3, which plots fjp against f\H for cores with varied magnetic 
fields and M = 10" 5 M yr -1 at 1 au and t = 10 4 yr. The 
direction of rotation changes with the sign of t)h, with the 
solid lines denoting negative rotation and the dotted lines 
positive. As the field increases the size of the area enclosed 
by the curve decreases, and the solutions tend towards the 
asymptote fjp = \t)h\ as both the diffusivities tend to zero. 
The ^p-intercept occurs at 



r\p = 



2b z i 



1+W1 + 



4m 2 



1/2 



(65) 



corresponding to the value of fjp needed to give the same 
magnetic field strength in KK02. The solid line intercepting 
the curves is fjp — | | a/ 2 + S 2 /S, which marks the bound- 
ary between the coefficient sets and in the limit of large |<5| 
also tends towards fjp = \fjn\- 

The second set of coefficients apply when the inequality 



> \8\ 



(66) 



is satisfied, so that B^s takes on the second value in Equation 
1581 B^s — —8B Z . The coefficient b Z 2 is given by the single 



„ rh 2 (l + 8 2 ) 6 m 6 
bz2 ~ 2/| & ^4/J-° 



where fa is 



h = fjp — t)h8\/2 + S 2 . 
This then gives the other coefficients as 

8bi 2 



n 



and 



(l + 8 2 )b 2 



-1 + 4/1 + 



4ra 2 



(l + <5 2 ) 2 6 4 2 



(67) 



(68) 



(69) 



(70) 



This solution was not explored in KK02, as without Hall 
diffusion the left term in Equation 1581 is zero. 

The second set of coefficients corresponds to the straight 
lines outside of the curves in Fig. 3, again with negative 
rotation drawn as solid lines and positive as dashed. The 
slope of the line (and the direction of rotation) are both 
dependant upon the sign of 8, i.e. the cap on B</, s ; in this 
solution the cap may have either sign, corresponding to the 
direction in which the field was azimuthally bending from 
the vertical before the cap was attained. Under a global 
reversal of the field the infalling matter spins in the opposite 
direction. 

The equations may also be rearranged to give the vari- 
ables as functions of the surface density and magnetic field: 



2GM 



M = -2nrY.V r t 
M = -2-KrY.Vr 



and J 



B(f, s B z 2 

2-K^Vr r ' 



(71) 

(72) 
(73) 

(74) 



where the field components are defined in Equations 
and the scale height of the disc is given by the equations 



#i = 
H a = 



hBlr 7/2 t 
V2GM 3 

(1 + 8 2 )B 2 z r 7 ' 2 t 
V2GM 3 



-i + Wi 



4c 2 M 2 



4c 2 , M 2 



{l + 8 2 ) 2 BjrH 2 



(75) 
(76) 



depending on which value of B^ a is adopted. 

This similarity solution exemplifies the magnetic brak- 
ing catastrophe that occurs in numerical simulations of star 
formation when the magnetic braking of a collapsing core is 
so strong that all angular momentum is removed from the 
gas and i t is impossible to form a rotationally-supported 
disc (e.g. iMellon fc L"Hl2008l , 120091 '). It is clear from Equa- 
tions !55b . l58l and l63l that when there is no Hall diffusion the 
magnetic braking causes B^g = J = 0, which would prevent 
Keplerian disc formation. 

The introduction of Hall diffusion to the similarity so- 
lution can cause additional twisting of the field lines and 
magnetic braking, or it can cause a reduction in magnetic 
braking by twisting the magnetic field lines in the opposite 
direction, in effect spinning up the collapse. The direction of 
the Hall diffusion depends upon the orientation of the field 
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with respect to the axis of rotation, and it is obvious that 
this directionality has an important effect on the magnetic 
braking catastrophe. The linear scaling of the angular mo- 
mentum with radius (Equation 1551) suggests that the point 
mass at the origin has no angular momentum, however, Hall 
diffusion will likely ensure that the angular momentum shall 
drop to a plateau value similar to that in the ambipolar 
diffusion-only collapse of KK02's fig. 10 in full simulations 
of such accretion flows. 



matching onto the power law solution described in f|3]at the 
inner boundary. Enforcing disc formation, Hall diffusion was 
shown to affect the mass of the protostar and accretion disc, 
as well as the size and structure of the disc. The free fall so- 
lution has yet to be adopted as a boundary condition in the 
full model, but its applicability to collapsing cores must be 
studied in future to show the importance of Hall diffusion 
in spinning up the flow, facilitating disc formation in order 
to solve the magnetic braking catastrophe. 



5 DISCUSSION 

The two similarity solutions represent both sides of the mag- 
netic braking catastrophe of star formation. In the first so- 
lution, the magnetic braking is limited and a rotationally- 
supported dis c simi lar to those in the simulations of 
iMachida et al.l (|201 lh forms, while in the second no disc 
forms as catastrophic magnetic braking remove s almost all of 
the an gular momentum as in the simulations of lMellon fc Lil 
(2009). The magnetic diffusion, in particular Hall diffusion, 
is clearly important in determining whether a disc forms, its 
size and rotational behaviour and the density profile within 
the disc. 

This can have wide implications on planet formation, 
for the size and surface density of the disc will then depend 
upon the product t]h(B • Q,), where the sign of r\H depends 
upon the abundances of charged grains and electron s rela- 
tive to the neutrals in the disc (jWardle fc Ndll999l ). This 
means that the outcome of planet formation in two other- 
wise identical discs can be vastly different if they are initially 
rotating in opposite directions. Our results in Sj3]show that 
one does not require a large amount of Hall diffusion for 
this asymmetry to be observable in the surface density of a 
protostellar disc. 

Our second solution shows a similar asymmetry caused 
by Hall diffusion, as all of the initial angular momentum 
has been removed from the collapse by magnetic braking 
and the only rotation of the flow is that induced by Hall dif- 
fusion. The direction of rotation corresponds directly with 
the sign of the Hall parameter, however the angular momen- 
tum remains dynamically unimportant. This solution ought 
to apply in simulations of accretion flows where the mag- 
netic braking parameter a and the cap on the azimuthal 
field 8 are large, as in the "strong braking" star formation 
similarity solution of KK02. 

The prescription for the magnetic braking used in this 
simple model is limited by our use of a cap placed on B^s 
to account for missing physics such as non-ax isymmetric ef- 
fects (e.g. the magnetorotational instability; ISano fc Stond 
2002). The azimuthal field takes the value of this cap in 
our Keplerian disc solution, and also in the free fall solu- 
tion when the Hall diffusivity is large, however the value 
5—1 adopted in our plots is that expected from a disc 
with a disc wind. A more complete description of the angu- 
lar momen tum transport would inclu de the calculation of a 
disc wind (jBlandford fc Pavndll982l ). as these are observed 
in simulations of protostellar and other accretion discs (e.g. 
IMachida et al.ll2007l : [Mellon fc Lillioolih. 

The similarity solutions of iBraiding fc Wardld (|2012l ) 
all demonstrated Keplerian disc formation in the inner re- 
gion of the protostellar collapse of a molecular cloud core, 



6 CONCLUSIONS 

We have presented two distinct power law similarity solu- 
tions to the MHD equations describing a flattened accre- 
tion flow with magnetic diffusion. The first of these repre- 
sented a rotationally-supported disc through which gas is 
slowly accreted while the magnetic field diffuses outwards 
against the flow. The second was a free fall collapse onto 
the star in which almost all of the angular momentum has 
been removed from the fluid and no rotationally-supported 
disc may form. These have been used to show that even a 
small amount of Hall diffusion can dramatically change the 
dynamics of gravitational collapse. 
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APPENDIX A: SELF-SIMILAR EQUATIONS 

Gravitational collapse occurs over many orders of magnitude 
in density, so that the point mass has negligible dimensions 
in comparison with the accretion flow. The self-similarity of 
the wave of infall, which propagates outwards at the speed of 
sound, arises because of the lack of characteristic time and 
length scales. Then the only dimensional quantities are the 
gravitational constant G, the isothermal sound speed c s , the 
local radius r and the instantaneous time t, so that, except 
for scaling factors, all quantities must be functions of the 
similarity variable x, defined as: 



These are augmented by the supplementary definitions: 



r 

C s t' 



(Al) 



We then define a set of nondimensional fluid variables that 
depend only upon x: 

E(r,t) = (^) *(x), g r (r,t) = (^) g(x), (A2) 
V r (r,t) = c 3 u(u), H(r,t) = c a th(x), (A3) 

V^r, t) = c.v(x), J(r, t) = c 2 tj(x), (A4) 

M(r, t) = m(x), M(r, t) = m(x), (A5) 

B ( r > *) = ( gTT^) b M, *^ = (fpr ) ( A6 ) 

In the thin disc approximation used here, the diffusivi- 
ties are assumed to be constant with height within the disc, 
and the Ohmic and ambipolar diffusion terms scale together 
to a zeroth-order approximation. They are combined into 
the Pedersen diffusivity, which scales within the disc as 



Vp =f]p- 



cH b 2 ^' 2 

T 3/2 



(A7) 



the self-similarity of the diffusivity depends upon the im- 
plicit relationship p oc pi, which holds t rue across a large 
range of densities dKamava fc Nish"ill2000l ) . We adopt a sim- 
ilar scaling for the Hall diffusion coefficient vh as a matter 
of pragmatism: 



Vh =f] H - 



c 2 t b 2 h^ 2 

7 3/2 ■ 



(A8) 



where vh is the nondimensional Hall diffusion parameter 
used to characterise the solutions, the ratio of the ambipolar 
to Hall terms becomes the most important factor in describ- 
ing the magnetic behaviour of the solutions. 

For convenience we define w = x — u and then use the 
similarity variables to write Equations [51 4131 as: 



dm 
dx 

(l- w 2 )±^L =g + 
a dx 

dj 



(A9) 



w 

X 



dx to 



^^h 2 + (b 2 s + b 2 s + a 2 )h - 2a — 0, 
x 

iji — xwb z + (fjub^^b + fjpb r ,sb z )xb z h 1/2 a~ 3 ^ 2 = 0, (A13) 

and — — = xb z . 
dx 



m = xwa, 
m — —xua, 

and g = -, 

x z 



while the field components are given by 



br.S — 



x 



mh>i 2 b 



a 3/2 

2arjph 1/2 ifjb z 



i— 

dx 



K 2_3/2 



(A15) 

(A16) 
(A17) 

(A18) 



(A19) 



a — a\x 

b z — b z ix~ q 



We find the similarity solutions by assuming that the 
fluid variables take the form of power law relations in x: 

(A20) 

(A21) 

(A22) 

where p, q, and r are real numbers. We substitute these into 
the self-similar equations, and by taking the limit as x — > 
we find the dominant terms and solve for the exponents and 
coefficients of the power law r elations. There are only two 
nontrivial solutions (derived in [Braiding 201l|): a Keplerian 
disc with p = 3/2, q — 5/4 and r = —1/2; and the near-free 
fall collapse with p = 1/2, q — 1 and r = — 1. 



Al Keplerian disc solution 

The first solution is the Keplerian disc presented in Sj3[ where 
the mass is supported against gravity by its angular momen- 
tum. The similarity solution in nondimensional form is: 

-3/2 



br,s — 

h - 



m 1/2 
-x b z , 



a — a±x 
j = ±V rhx, 



3*' 
h lX s ' 2 , 



b z = 
b<t>,3 = 



in 



3/4 



/2|« 

-Sb z 



z.T 



-5/4 



(A23) 
(A24) 

(A25) 

(A26) 
(A27) 



where the constants are defined as in ij3] In the ambipolar 
diffusion limit these reduce to equations 51-57 of KK02. 

A2 Free fall solution 

The second solution is the set of relations that describe a free 
fall onto the central mass dominated by magnetic diffusion: 



(A10) 


m = 


m, 


a 


j rh 
= \2x-' 




(A28) 


(All) 


u = 


[2rii 

IT' 


3 


= hfix, 




(A29) 


(A12) 




b z x, 


b z 


= b z ]_,2X~ 


-l 


(A30) 


(A13) 


br,s = 


bz, 




— — min 




, (A31) 












.VP 


(A14) 


h = 


l, 3/2 








(A32) 
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Given that there are two possible values for b</, jS) depending 
on the chosen values of the parameters fjn,A and S, there are 
two sets of coefficients for j, b z and h. These are defined in 
fj4] and referenced by the subscripts f and 2 depending on 
which term in Equation IA3fb they are associated with. 

In the ambipolar diffusion limit this solution reduces to 
equations 66-71 of KK02, where j = 6^ iS « 0. 

This paper has been typeset from a TgX/ I^TjjX file prepared 
by the author. 



The Hall effect in accretion flows 9 



© 2012 RAS, MNRAS 000, [TH] 



